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Abstract 

In this work we consider a two-dimensional piecewise smooth system, 
defined in two domains separated by the switching manifold x — 0. We 
assume that there exists a piecewise-defined continuous Hamiltonian that 
is a first integral of the system. We also suppose that the system possesses 
an invisible fold-fold at the origin and two heteroclinic orbits connecting 
two hyperbolic critical points on either side of x = 0. Finally, we assume 
that the region closed by these heteroclinic connections is fully covered 
by periodic orbits surrounding the origin, whose periods monotonically 
increase as they approach the heteroclinic connection. 
When considering a non-autonomous (T-periodic) Hamiltonian perturba- 
tion of amplitude e, using an impact map, we rigorously prove that, for 
every n and m relatively prime and e > small enough, there exists a nT- 
periodic orbit impacting 2m times with the switching manifold at every 
period if a modified subharmonic Melnikov function possesses a simple 
zero. We also prove that, if the orbits are discontinuous when they cross 
x = 0, then all these orbits exist if the relative size of e > with respect 
to the magnitude of this jump is large enough. 

We also obtain similar conditions for the splitting of the heteroclinic con- 
nections. 

Keywords: subharmonic orbits, heteroclinic connections, non-smooth im- 
pact systems, Melnikov method. 



1 Introduction 

The Melnikov method provides tools to determine the persistence of periodic 
orbits and homoclinic/heteroclinic connections for planar regular systems under 
non-autonomous periodic perturbations GH83 . This persistence is guaranteed 
by the existence of simple zeros of a certain function, the subharmonic Melnikov 
function and the Melnikov function, respectively. In this work we extend these 
classical results to a class of piecewise smooth differential equations generalizing 
a mechanical impact model. In such systems, the perturbation typically mod- 
els an external forcing and, hence, affects a second order differential equation. 
However, in this work, we allow for a general periodic Hamiltonian perturbation, 
potentially influencing both velocity and acceleration. Note that no symmetry 
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is assumed in either the perturbed or unperturbed system. 

The unperturbed system is defined in two domains separated by a switching 
manifold where the impacts occur, and possesses one hyperbolic critical point 
on either side of it. We distinguish between two situations regarding the un- 
perturbed system. In the first one, which we call conservative, two heteroclinic 
trajectories connect both hyperbolic points, and surround a region completely 
covered by periodic orbits including the origin. In the second one, we introduce 
an energy dissipation at the impacts, which is modeled by an algebraic condition 
that forces the solutions to undergo a discontinuity every time they cross the 
switching manifold. Then, the origin becomes a global attractor and none of 
these objects can exist for the unperturbed system. 

For a smooth system, the classical Melnikov method considers fixed (or pe- 
riodic) points of the time T stroboscopic map, where T is the period of the 
perturbation. However, for our class of system, such a map becomes unwieldy 
because one has to check the number of times that the flow crosses the switching 
manifold, which is a priori unknown and can even be arbitrarily large. Instead, 
using the switching manifold as a Poincare section and adding time as variable, 
we consider the first return Poincare map, the so-called impact map. This map is 
smooth and hence we can use classical perturbation theory to rigorously prove 
sufficient conditions for the existence of periodic orbits. In the conservative 
case, these conditions turn out to be same ones given by the classical Melnikov 
method, so extending it to a class of piecewise smooth systems (Theorem!]}. In 
addition, we rigorously prove that the simple zeros of the subharmonic Melnikov 
function can guarantee the existence of periodic orbits when the trajectories are 
forced to be discontinuous due to the loss of energy at impact (Theorem [5]) . 
In addition, the impact map could also be used to prove the existence of invari- 
ant KAM tori in the system. After writing the system in action-angle variables, 
these ideas were applied in |KKY97| to a different system to prove the existence 
of such tori. 

To prove the existence of heteroclinic connections for the perturbed case, it 
is sufficient to look for the intersection with the switching manifold of the stable 
and unstable manifolds |BK91j . |Hog92| . In this way, we rigorously extend the 
classical Melnikov method for heteroclinic connections to a class of piecewise 
smooth systems. When the loss of energy is considered, we prove that the ze- 
ros of the Melnikov function guarantee the existence of transversal heteroclinic 
intersections. Both results are given in Theorem [3] 

This paper is organized as follows. In Sj5]we describe the class of system that 
we consider, state some notation and introduce tools needed for this work. In Sj21 
we prove the existence of periodic orbits distinguishing between the conservative 
and dissipative cases. SjHis devoted to heteroclinic connections. Finally, in we 
use the example of the rocking block to illustrate the results obtained regarding 
the periodic orbits, and compare with the work of [Hog89 . 
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2 System description 

2.1 General system definition 

We divide the plane into two sets, 

5+ = {(.T,y)eM 2 |x>0} 
S~ = {0,y) G M 2 \x < 0} 

separated by the switching manifold 

S = £+ UE - U (0,0) (2.1) 

where 

S+ = {(x,y)e]R 2 |2;>0} 

sr = eR 2 \x < o|. 

We consider the piecewise smooth system 

f X+{x,y)+eX+{x,y,t) ii (x,y) G S+ 



a; 



2// I # (x, y) + eX x (x, y, t) if (x,y) G S~ 



(2.2) 



We assume G C°°(IR 2 ) and Aff^y, f) G C°°(E 3 ), although this can be 
relaxed to less regularity in S ± and x R, respectively. 

System (|2.2p is a Hamiltonian system associated with a C° piecewise smooth 
Hamiltonian of the form 

H e (x,y,t)=H (x,y)+eH 1 (x,y,t). (2.3) 
The unperturbed C°(1R 2 ) Hamiltonian Ho is a classical Hamiltonian given by 

2 

I H+ (x, y) := + V + (x) if (x, y) G 5 + U S 



ff (*,y) :=^-+V(:r) := 



H -(a;,y):=^-+V-(a:) if(a;,y)eS- 



(2.4) 

with V ± G C*°°(IR 2 ) satisfying V+(0) = V~(0). 

Similarly, the non-autonomous T-periodic C°(R 3 ) perturbation, eH\, is given 

by 

H+(x,y,t) if (x,y) G5+UE+ 
Hi(x,y,t) if (x, y) G S" - 

fulfilling Ht(0,y,t) = H^(0,y,t) V(y,t) G M 2 . 

Then, the relation between (12. 2j) and (|2.3[) is given by 



Hi(x,y,t) 



X++eX+ = JV(H++eH+) 
X - + = JV(ff _ + eHT), 
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where J is the usual symplectic matrix 



J = 



( 







I 1 
1 



) 



We assume that the phase portrait of the unperturbed system (|2.2I) (e = 0) is 
topologically equivalent to the one shown in Fig. [1] which we make precise in 
the following hypotheses. 

C.l There exist two hyperbolic critical points z + = (x + ,y + ) G S + and z~ = 
(x~,y~) G S~ of saddle type belonging to the energy level 



C.2 The origin is an invisible fold-fold of centre type [GST11] . such that Hq(0, 0) = 



C.3 There exist two heteroclinic orbits given by W u (z ) = W s (z + ) and 
W u {z + ) — W s (z~) surrounding the origin and contained in the energy 
level (I2~lj) . 

C.4 The region between both heteroclinic orbits is fully covered by periodic 
orbits surrounding the origin given by 



with < c < ci, and A c intersects £ transversally exactly twice. 

C.5 The period of A c is a regular function of c with strictly positive derivative 
for < c < ci . 

Note that, as the unperturbed Hamiltonian Hq is C°° in S + and S~ , the fact 
that the heteroclinic orbits are in the energy level Ho(x,y) = c\ follows auto- 
matically from hypothesis C.l. However, we include it explicitly for clarity. 

We wish to determine which of these objects and characteristics persist and 
which are destroyed when the small non-autonomous T-periodic perturbation 
eHi is considered. Of interest is the splitting of the separatrices and the persis- 
tence of periodic orbits. In the smooth case, these answers are given completely 
by the classical Melnikov method [GH83] . Hence, it is natural to check whether 
these classical tools are still valid for the piecewise smooth system presented 
above and if any changes to the method are necessary. 

Another interesting question that can be addressed with a similar approach 
is the existence of 2-dimensional invariant tori of system (|2.2p (see [KKY971 




(2.6) 



0. 




(2.7) 



IKunOOp . 
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Figure 1: Phase portrait for the unperturbed system (|2.2[) . 
2.2 Poincare impact map 

To study system (|2.2j) we will proceed as in |Hog89| using the Poincare impact 
map. We consider the extended phase space l 2 xl adding time as a system 
variable and equation t = 1 to Eq. (|2.2|) . As the perturbation is periodic, this 
time variable is usually defined in T = R/T; however, it will be more useful for us 
to consider R instead. We want to study the motion in the region surrounded by 
the heteroclinic orbit, so we consider in this extended phase-space the Poincare 
section 

S+ = {(0,y,t) G R 2 x R|0 < y < V^ci} ■ (2.8) 

To simplify the notation, as the first coordinate in S + is always 0, we will omit 
its repetition whenever this does not lead to confusion. The domain of the 
Poincare map is not S + but a suitable open set U , that depends on e and, for 
e = 0, docs not contain the heteroclinic connection. 

We now define the Poincare impact map 

P £ :U <zt+ — > E+, 
as follows (see Fig. First, using the section 

XT = {(0,y,t) e R 2 x R| - y/2c[ < y < 0} , (2.9) 
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with (0, yo, to) e U + C E + , we define the map 

P+ : U + C £+ — > IT, 

as 

P+(W),*o) = (n»(^ + (*i;to,0,|«),e)),ti) (2.10) 

where <p + (t; to , x , y , e) is the flow associated with system (|2.2[) restricted to 
and t\ > to is the smallest value of t satisfying the condition 

n*fa + (ti;to,o,i/o,e))=o. (2.11) 

Similarly, we consider 

P~ : U~ C XT — ► E+ 
for (Q,j/i,*i) G C E~ defined by 

iT(l/i.*i)=(n v (0~(t2;ti,O,v 1 ,e)),t2) (2.12) 

where (j}~(t;ti,x,y,e) is the flow associated with (|2.2I) restricted to S* - , and 
*2 > t\ is the smallest value of t satisfying the condition 

n*(0~(t2;ti,O,ift,e))=Q. (2-13) 
Then the Poincare impact map is defined as the composition 
P £ : U C E+ -> E + 

(j/0,to) ' ^ ProP/feo.io) 

Notice that, as assumed in C.4, for the unperturbed flow all initial conditions in 
E + lead to periodic orbits surrounding the origin. Hence, we can give a closed 
expression for Pq, the Poincare impact map when e = 0. Let 

a ± (±y)=±2[ ■ 1 dx, h = H (0,±y) = — (2.15) 

1 ' J J2{h-V ± {x)) K y ' 2 y ' 

o 

be the time needed by an orbit of the unperturbed system with initial condition 
(0,±y) G to reach E T . In the unperturbed case, the orbit with initial 
condition (0, y) G E + has period 

a(y)=a + (y) + a-(-y). (2.16) 

Then the Poincare impact map when e = is defined in the whole E + , and can 
be written as 

Po(yo,to) = (yo,to + a(y )). (2.17) 

Thus, if £ is small enough, the perturbed trajectories starting at E + cross E + 
again. The Poincare impact map is well defined, and is as smooth as the flow 
restricted to S + and S~ . 

Note that in the symmetric case, V + (x) = V~(— x), a + (y) = a~(—y) is half 
the period of the unperturbed periodic orbit with initial condition (0,y) G E + . 
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(yo,to) 
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Figure 2: Poincare impact map ()2.14|) represented schematically. 
2.3 Coefficient of restitution 

As the name of the previous map suggests, it is typically used to deal with 
systems with impacts, as is the case of the mechanical example of section [5j 
In order to include the loss of energy at the impact, one considers a coefficient 
of restitution, r £ (0,1], that reduces the velocity, y, at every impact. More 
precisely, if a trajectory crosses E transversally at some point (0,ys) at t = ts, 
then the state is replaced by (0,rj/s) at a later time tA to proceed with the 
evolution of the system. In other words, the system slides along S from (0, ys) 
to (0,t-?/b) during time tA — is and 



For the rest of this article we will assume that the loss of energy is produced 
instantaneously and hence tA = is- Thus, there is no sliding along E and the 
trajectory jumps from (0, ys) to (0, rys)- 

Clearly, when such a condition is introduced to a system of the type f|2 . 2[) . 
the unperturbed system (e = 0) is no longer conservative, the origin becomes a 
global attractor and none of the conditions C.1-C.5 hold. In particular, the or- 
bits with initial conditions on the unstable manifolds W u (z~) and W u (z + ) tend 
to the origin and can not intersect the stable manifolds W s (z + ) and W 3 (z~), 
respectively (see Fig. [3]). 

Although periodic orbits surrounding the origin are not possible for the un- 
perturbed case if r < 1, they may exist if e > 0. However, roughly speaking, 
as these orbits will have to overcome the loss of energy, the magnitude of the 
forcing will not be allowed to be arbitrarily small. We will make a precise state- 
ment of this fact in §3.21 fsee also |Hog89 ). 



v{tA) = ry(t B ). 



(2.18) 
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To study the existence of periodic orbits we will use again the impact map, 
which can also be defined for r < 1 as (see Fig. 0} 

PeAyo,to) ■= Rr o P~ o Br o P+(t/ ,*o) (2-19) 

where 

Br(yo,to) = (ry ,t ). 

Note that P e>r is as smooth as the flow restricted to S^, since it is the compo- 
sition of smooth maps. 

Using Eqs. (|2.15[) and (|2.16[) . the impact map, P £)r , for e = and r < 1 can 
be written as 

pQ.Ayo-.to) = (r 2 y ,t + a + (y ) + a-(-ry )) . (2.20) 
Note that, for any e > 0, 

P £ ,i(yo,*o) = P e (yo»<o)- 
2.4 Some formal definitions and notation 

Up to now, we have considered separately the solutions of system (|2.2j) in S + 
and S^until they reach the switching manifold E. Given an initial condition 
(xo,yo,to), one can extend the definition of a solution, (j>(t;to 7 xo,yo,£,r), of 
system (|2.2p . (|2.18|) for all t > to by properly concatenating <fi + or <fi~ whenever 
the flow crosses E transversally. Depending on the sign of Xq, one applies either 
(j) + (t;to,xo,yo,e) or (f>~ (i; to, xq, yo, e) until the trajectory reaches E, and then 
one applies (|2.18j) . If xq = 0, one proceeds similarly depending on the sign of 
yo- This is because x — y + 0(e) is always an equation of the flow and the orbits 
twist clockwise. 

In this work, we will mainly use solutions with initial conditions (0, yo, ^o) € 
E + . In that case, we define the sequence of impacts (0, yl r ,t l er ) (see Fig. [5]), if 
they exist, as 

f .,i a ,\RrO p-{yi-\t-i) if y\f < 

Ufe,r> t e ,r) ~\ RrQ P+(^-l,4-l) if y-1 > ' ^- ll > 

with (j/°r!*°r) = (Vo,to) and P± defined in (|2~T0| and (|2T2]> , Notice that the 
sequence (|2.21|) will be finite if the flow reaches E a finite number of times only. 
For the unperturbed case, for any point (0, yo,to) £ E + , the sequence (|2.21|) 
becomes 

,i a \ ._ / {r l yoA~r + (-r i_1 2/o)) if i > 2 even 
^o,r:VJ-| (-rVc^ + a+^yo)) if i > 1 odd ' [ ^ 2) 

where ar^ are defined in Eq. (12 . 15[) . 

Once the impacts (y\ r , t\ r ) are defined, the solution of the non-autonomous sys- 
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(y ,t ) = (y°e,r,tl) y 










X I * 







Figure 5: Sequence of impacts for r < 1 and e > 0. 

tem 

(|2.2[) . (|2.18j) with initial condition (0, yo,to) € S + is given as 

f 0+(i;t* ,0,i# r ,e) if^ r <t<^+i 
^(i;io,Q,j/Q,e,r) := < , i > 0. 

(2.23) 

Note that in the case when the number of impacts is finite, we take the last 
interval of time to be infinitely long. 

In the rest of the paper we will generally distinguish between the conserva- 
tive (r = 1) and dissipative (r < 1) cases. We will omit the parameter r in the 
flow </> whenever we refer to r = 1. 

Note that we have only defined the solution of the system for an initial condition 
(0, yo,to) € S + . Given (0,y O7 t o ) € £~, one defines similarly this solution by 
just properly shifting the subscripts of t\ in (I2.23|) . In addition, it is possible to 
extend precisely this definition to an arbitrarily initial condition (2^, yo, to). 

As is usual when dealing with Hamiltonian systems, we will use the unper- 
turbed Hamiltonian to measure the distance between states. In addition, as we 
are dealing with a perturbation problem, we will frequently use expansions in 
powers of e. In this case, the integral of the Poisson brackets of the Hamiltoni- 
ans Hi and Hq typically provides a compact expression for the linear terms in e. 
Given m > 1, (0, yo, to) € and its impact sequence (0, y\ r , t\ r ), < i < 2m, 
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for the non-smooth system (|2.2p . (|2.18p when r < 1, we introduce 



{H , Hi) ((/> (t; t , 0, y , e, r) , t) dt 



m-l / r t 2i+1 



i=0 \ 

t 2i + 2 



: = E / {H+,H+}(cj ) +(t;t 2 E : r ,0,y^ r ,e),t)dt (2.24) 



i' 2 



+ / {ff - )J ffr}(^-(t;^+ 1 ,0,t/^+ 1 , e ),t) 

/ , 2i+ 1 



rff 



where {Q (x, y) , R (x, y)} — §§§^ — T^ff is the usual Poisson bracket of the 
Hamiltonians Q an d R- 

The next Lemma provides an expression for Ho (4> (t 2 ™; to, 0, yo, £, r) ) which 
we will use below. 

Lemma 1. Let m > 1 and (0, y , t ) € £ + , and Zet (0, r , t\ r ), i = 0,. . . , 2m, 
be its associated impact sequence as defined in &2.21\) . Then, 



Ho(0,y 2 ™)-H (0,yo) = r 2 



e I {H o ,H 1 }((t)(t;to,0,yo,e,r),t)dt 



it, 



E \Ho(0,yi r )-H 0,^1 



(2.25) 



Proof. The proof of this Lemma comes from a straightforward application of the 
fundamental theorem of calculus to the smooth functions H^ (0 ± (i; to, 0, ±y , e)), 
using the fact that 

H Q (0, y 2 ™) = Ho {rcj> (t 2 ™;t , 0, yo, e, r)) 

= r 2 H (cj)(t 2 ™;to,0,yo,£,r)) , 

taking into account the intermediate gaps induced by the impact condition (12.181) 
and using the fact that 



jH± (i; t*, x* , y* , e)) = e {iJ ± , Hf } (t; f, , y*, e)) 



for any (x*,y*) € 5 U £ and f > t* such that (i; t* ,x* ,y* ,e) € S*. □ 

The following Lemma gives us an expression for the expansion in powers of 
£ of # (Q, - #o (0, y ) which we will use in [j3J 



Lemma 2. Lei m > 1 and (0, yo, to) € S + , and fei (0, r , i| r ), i — 0, . . . , 2m, 
be its associated impact sequence as defined in (2.21)) . Then, if £ ~ 0, the Taylor 
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expansion of expression \2.25\) becomes 



H (0, y 2 E m r ) - H (0, y ) = ^-(r im - 1) + eG m (y ,t ) 
+ 0(e 2 ) + 0(e (r-1)) 



wher 



G m (yo,io) = 



{# o ,/ii}(0 (t;0,0 )W) ,0),t + fo)dt 



(2.26) 



(2.27) 



and a(yo) is given in \2.16\) . 

Proof. The independent term of the expansion is found by noting that, if e = 0, 

from expression (|2.22l) . one has H (0, j/q r ) = H a (0 : -^f-). Hence all the terms 
in the sum of Eq. (|2.25[) cancel each other except for the first and the last one. 

2 

This, in combination with the fact that Ho(0,y) = H^-, gives the first term in 
Eq. (|2.26p . For the linear term in e, one first obtains 



{iTo,£ri}(0(t;*o,O > |/o,O,r),i)d* 



2m- 1 



(- 2 - 1 ) E (i(tf°((<^y))i £ =o 



Then, by applying m times Eq. (|2.20[) . one has that t*™ = t a + ma(y ) + 
0(e). Thus, by expanding this for r near 1 and e near and noting that the 
unperturbed flow is autonomous and hence <f>(t; to, 0, yo, 0) = <f>{t — to] 0, 0, yo, 0), 
one gets expression (I2.27[) . □ 

Remark 1. If in Eq. we take a(yo) — then we recover the clas- 

sical Melnikov function for the subharmonic orbits 'GH831 with the modified 
integral [E£$ . 



3 Existence of subharmonic orbits 

3.1 Conservative case, r = 1: Melnikov method for sub- 
harmonic orbits 

Let us consider system fl2.2[) neglecting the loss of energy at impact (r = 1 in 
Eq. (|2.18|l ). According to C.1-C.5, for e = 0, this system possesses a continuum 
of periodic orbits, A c in Eq. (|2.7[) , surrounding the origin. Our main goal in this 
section is to investigate the persistence of these orbits when the (periodic) non- 
autonomous perturbation is considered (e > 0). The classical Melnikov method 
for subharmonic orbits, which here, in principle, does not apply, provides suffi- 
cient conditions for the persistence of periodic orbits for a smooth system with 
an equivalent, smooth, unperturbed phase portrait. 
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The period of the orbits A c tends to infinity as they approach the heteroclinic 
orbit. More precisely, if q c (t) is the periodic orbit satisfying q c (0) — (0,yo) with 
-ffo(0,yo) = c, its period a{yo) tends to infinity as c — > c\ (see formula (|2.16p ). 
As we are interested in finding periodic orbits for < e <C 1, we will use the 
unperturbed periodic solutions as e-close approximations to them. In general, 
such perturbation results are valid only for finite time and therefore, from now 
on, we will restrict ourselves to a set of the form 

S c + = {(0,y,<)eS+|0<y<c}, (3.1) 

for a fixed c satisfying < c < \f2c[. Note that, if (0, yo^o) G then a{yo) 
is uniformly bounded (a(yo) < <*(£)). However, following GH83 , it is also 
possible to extend the method for all the periodic orbits up to the heteroclinic 
connection. 

To look for periodic orbits we will use the impact map defined in (|2.14|) . In 
terms of this map, a point in U C S + will lead to a periodic orbit of period nT 
if it is a solution of the equation 

P™(y ,t Q ) = (y ,t +nT), (3.2) 

for some m. We take m to be the smallest integer such that (|3.2[) is satisfied. 
In that case, 4> (t] to, 0, yo, s) will be a periodic orbit of period nT, which crosses 
the switching manifold £ exactly 2m times. We will call this an (n, m)-periodic 
orbit. Then for (n, m)-periodic orbits with e > we have the following result 
analogous to the smooth case 

Theorem 1. Consider a system as defined in h2.2\) satisfying C.1-C.5, and let 
ot{yo) be the function defined in H2.15\) - fr2.16\) . Assume that the point (0, yo, to) G 
T,^ satisfies 

H.l a(yo) — ^-T with n,m S Z relatively prime 
H.2 to G [0,T"] is a simple zero of 

1-nT 

M n > m {to)= {Ho,Hx}{q c {t),t + h)dt, c = H (0,y ), (3.3) 
Jo 

where q c (t) = 4>(t;0, 0,y~o) is the periodic orbit such that ct(y~o) = 7^"- 

Then, there exists £0 such that, for every < e < £0, one can find z/q and ig 

such that (f>(t;tQ,0,yQ,e) is an (n,m) -periodic orbit. 

Proof. The proof of the result comes from a straightforward application of the 
implicit function theorem to equation (I3.2[) . Let us fix n and m relatively prime. 
We replace equation (|3.2I) by 

/ H (0,IL yo (Pp(yo,to))) \_f H o (0,yo) \_ ( A 

1, II to (P e "%o,to)) ) {t o + nTj-{0j- {6A) 
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That is, we use the Hamitlonian Hq to measure the distance between the points 

(0,n yo (P™(y ,t ))) and (0,y ). 

Using the second equation in (|3.4I) we have 

U yo (P™ (y , t )) = U y (0(t o + nT; i„, 0, y , e)) 
= n x (0(t o + nT;t Q ,O,y o ,e)). 



This allows us to rewrite Eq. (|3.4I) as 



H (4>(t + nT;to,0,yo,s))-H (0,yo) \ ( 
n to (P E ro (yo,to))-nT-t / lo 



(3.5) 



We expand Eq. (I3.5[) in powers of e. Using Eq. (|2.17p . the second component 
of (|3.5[) becomes 

n t0 (P £ m (y , t )) -h-nT = ma{y ) - nT + 0(e) = 0, (3.6) 
where a(yo) is the period of the periodic orbit q e (t), c — Ho(0,yo), given in 

Eq. ipnsji . 

On the other hand, using Lemma [5] and noting that 

ll yo (Pr(yo,to)) = y 2 J, 

the first equation in ((3 . 5|) can be written as 
Ho(0,Tl yo (P™(yo,to))) - H o (0,yo) 

rma(yo) 

= e / {Ho, Hi) {cf)(t;0,0,y ,a),t + t )dt + O{e 2 ) 

Jo 

= eG m (yo,to) + 0(e 2 ). 
where G m (yo,to) is given in (|2.27|l . Hence, Eq. (I3.5[) finally becomes 

F n , m (y 0) t ,e) := ( ^J^l^+o^) ) = ( ) ' (3J) 

where the order in e of the first component has been reduced and, thus, the 
implicit function theorem can be applied to Eq. (|3.7p . Therefore, one needs 

1. F ntm (y ,i o ,0) = (0,0) T 

2. det(D yo j F(yo, to, 0)) ^ 0, where D VOi t = D is the Jacobian with respect 
to the variables yo and to- 

The first condition is satisfied by noting in Eq. (|3.7[) that yo has to be such that 
a (yo) — 7^" and to a zero of the subharmonic Melnikov function 

MlT 

M"< m (t ) :=G m {yoM) = / {H , H x } (q c (t),t + t )dt, 
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where q c (t), c = H o (0, y ), is the unperturbed periodic orbit of period ^ such 
that q c (0) = (0,j/o)j an d therefore q c (t) — <j> (t; 0, 0, yo, 0). 
In addition, for e — 0, DF n ^ m is given by 



V ™" (2/o ) 



By C.5, a'(yo) ^ 0, and the second condition is satisfied if to is a simple zero 
of the subharmonic Melnikov function, M n,m (to), which completes hypothesis 
H.2. 

Finally, applying the implicit function theorem to (I3.7[) at (y , to,s) = (yo, to,0), 
there exists £o > such that, if < e < So, then there exist unique yo(s) and 
^o(e) solutions of the equation (|3.4[) . which have the form 

2/o = So + O(e) 
t* = t + 0(e). 

Hence, the orbit <fi (t; t^, 0, y^, e) is an (n, m)-periodic orbit, as it has period nT 
and impacts 2m times with the switching manifold S in every period. □ 

Remark 2. The upper bound Eq given in the theorem depends on n and m. 
However, for every fixed m, it is possible to obtain So(m), such that for e < 
eo(m), we can apply the theorem for all n such that a -1 (^-) € Eg . TTiiis is 
because the approximation of the perturbed flow by the unperturbed periodic orbit 
is performed m times beyond the period of the unperturbed periodic orbit. 

Remark 3. The proof of the result provides us with a constructive method to 
find the initial condition for nT -periodic orbits for e > 0. 

1. Given n and m, find yo such that a(y~o) = using Eq. 12.16\) . 

2. Find to such that M n,m (to) has a simple zero at to = to- 

3. Use (yo,io) o,s seed to solve Eq. \3.$ numerically. 

Lemma 3. The subharmonic Melnikov function US. 3\) is either identically zero 
or generically possesses at least one simple zero. 

Proof. The proof comes from the fact that M n ' m (to) has average 

1 f T 

1 Jo 



equal to zero. 



1 f T r T 

<M n < m (to)> = - / {Ho,H 1 }(q c (t),t + to)dtdt 



/ / {Ho,#l}(9c(<),t 

Jo Jo 

/ {Ho,<H x >}(q c (t))dt. 
Jo 



to) dtodt 
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Recalling that a(y ) = ^- (see ^J^-^J^) and letting 

q+(t) = <t>+(t; 0,0,2/0,0) 

£T(*) = ^(i;" + (yo), 0,-yo,0), 

< Af"' m (t ) > can be written as 

f-a + (j/o) 



/ ° {H+, < H+ >} (g+(t)) dt + f m {Ho, < >} (?j(t)) dt ) 

JO Ja+(y ) ) 

(< Jf+ > ( ? +(t))) + / I (< H f > (g"(t))) dt 



-m 



dt ' ' ' Ja+(y ) dt 



= -ml < H+ > (qt(a+(y Q ))) - < H+ > (q+(0)) 
+ < #r > (g-(^)) - < H x - > («+(y ))) ) = 0. 



□ 



Note that, if M"' m (t ) = then a second order analysis is required to study 
the existence of periodic orbits. 



3.2 Dissipative case, r < 1 

We now focus on the situation when the coefficient of restitution r introduced in 
£|2.3l is considered. As already mentioned, for e = the origin is a global attrac- 
tor and hence none of the periodic orbits studied in the previous section exists 
if the amplitude of the perturbation is small enough. However, as was shown 
in |Hog89| for the rocking block model, for e large enough an infinite number 
periodic orbits surrounding the origin can exist. This was studied analytically 
and numerically for the rocking block model under symmetry assumptions for 
the particular case m = 1. Here, our goal is to relate the periodic orbits exist- 
ing for the dissipative case to those which exist for r = 1 in the general system 
(|2.2I) . (|2.18[) . As will be shown below, all the periodic orbits given by Theorem [T] 
can also exist for the dissipative case, when r < 1 is small enough compared 
with e > 0. In other words, we generalise in this section the result presented 
for the conservative case. 



As in i j3.ll in order to obtain the initial conditions of a (n, m)-periodic orbit 
for r < 1, one has to solve the equation 

P™(y ,to) = (y ,t + nT), (3.8) 

where P £yr , is defined in Eq. (|2.19[) . The next result states that, under certain 
conditions regarding r and e, Eq. (|3.8I) can be solved. 
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Theorem 2. Consider system i2.2])J2.18\) . Let (0, yo,t ) e S + be such that 
a (Vo) = wit/l n and m relatively prime, and to a simple zero of the sub- 
harmonic Melnikov function i3.3\) . There exists p such that, given e, f > 
satisfying < ? < p, t/iere exists 6 such that, if e = eS and r = 1 — f <5, t/ien 
V0 < (5 < So there exists (j/oi^o) which is a solution of Eq. H3.8\) . Moreover, 
Vo — + 0(6), t^ — to + 0(6) + 0(f/e) and the solution (yo\^o) tends to the 
one given in Theorem^ as r — > 1 — . 

Proof. As in the conservative case, we use the unperturbed Hamiltonian to 
measure the distance between points in E. Then, Eq. (|3.8j) can be rewritten as 



H (Q,n yo (Pp(yo,t ))) \ / Ho(0,y ) \ = ( 
n tD (p™(y ,i )) J { t a +nT 







(3.9) 



As in Theorem [TJ we proceed by expanding this equation in powers of e and 
r - 1 using ([2726]) and (f2T27|) obtaining 



/ *f(r 4 ™ - 1) + eG m (yo,io) + 0(e 2 ) + 0(e(r - 1)) 



2 

m — 1 rn — 1 



£ « + (r 2 Vo) + E a ~ (-^ +1 ?/o) + 0(e) - nT j V 

V i=0 i—0 







(3.10) 



Note that, for r = 1, Eq. ([3"7T0]) becomes Eq. (|3"3j) . 

We are interested in studying Eq. (|3.10l) when 1 — r and e are both small. 
Therefore, for e > and f > we set 

s = e<5, r = 1 — r6, (3-H) 

where (5 > is a small parameter. Then Eq. p.!0[) becomes 

F n ,m(yo,t ,6) : = 

-2mry 2 + sG m (y , t ) + 0(6) \ _f0\ ^ 



ma(y ) + 0(6) - nT J \ 

We now need to apply the implicit function theorem to (|3 . 1 2[) . 

The first step is to solve Eq. (|3.12j) for 6 = 0. The second equation gives 

a (yo) — 7?" ; as m Theorem [TJ To solve the first equation, we define 

P> m (to) = -2mff + eM n > m (t ), (3.13) 

and to will be given by a zero of f n,m (to)- Assume now that to is a simple zero 
of M n ' m (to). As M n ' m (to) is a smooth periodic function, it possesses at least 
one local maximum. Let tu be the closest value to to where M n ' m (to) possesses 
a local maximum, and assume (M n ' m )' (to) ^ for all to between to and tM- 
If (M n ' m )' (t ) vanishes between t and tjv/, we then take tu to be the closest 
value to t such that (M n ' m )' (t ) = to ensure that (M n,m )' (t ) ^ between 
to and t M . We then define p := M "" 2% m) ■ Then, if 

zmyQ 

< ~ < p, (3.14) 

£ 
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there exists to |-close to to where f n ' m (to) has a simple zero. Since a'(yo) > 0, 
a similar calculation to the one in Theorem Q] shows that 

dot (pP yoM (£ ,*o,0)) ^0, 

and hence we can apply the implicit function theorem near (yo, to, 5) = (yo, to, 0) 
to show that there exists So such that, if < 5 < So, then there exists 

(VoXo) = (vo,%) + 0(S) = (y ,i ) + O(S) + O(r/g) 

which is a solution of Eq. (|3.8I) . 

This solution tends to the one given by Theorem [1] when r — > + . This is a 
natural consequence of that fact that Eq. (|3.8|) tends to the Equation (|3.2j) as 
r4l". □ 

Remark 4. In order to determine p, we have imposed tu to be the local maxi- 
mum of the Melnikov function closest to its simple zero, to- Instead, one could 
also use the absolute maximum so increasing the range given in Eq. I[3.14\ ). 
However, in this case, the values where (M™'™ 1 )' (to) = have to be avoided to 
ensure that the desired zero of f n ' m (to) is simple. 

Remark 5. Arguing as in Remark^ for every m fixed, the constant So(m) can 
be taken such that if S < So(m), there exist periodic orbits for all n such that 



4 Intersection of the separatrices 

We now focus our attention on the invariant manifolds of the saddle points of 
system (|2.2I) . (|2.18I) when e > 0. As explained in ^ for e = 0, there exist two 
heteroclinic orbits connecting the critical points 2F if r = 1 (see Fig. [TJ whereas, 
if r < 1, the unstable manifolds W u (z ± ) spiral discontinuously from z ± to the 
origin and W s (z ± ) becomes unbounded (see Fig. [3]). As we will show, in both 
cases, heteroclinic orbits may exist for the perturbed system. 

For a smooth system with Hamiltonian Ko(x,y) + eKi(x,y,t), the persis- 
tence of homoclinic or heteroclinic connections is achieved by the well known 
Melnikov method which states that if the Melnikov function 

/+00 
{K ,Ki}(4>(t;t ,z ,0),t + t )dt, 
-00 

with zo — (xo,yo) £ W u (z~) = W s (z + ), has a simple zero, then the stable and 
unstable manifolds intersect for e > small enough (see [GH83 ). 
In this section we will modify the classical Melnikov method and we will rigor- 
ously prove that it is still valid for a piecewise smooth system of the form (|2.2[) , 
even if r < 1 . 
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Figure 6: Section of the unperturbed and perturbed invariant manifolds for 
t = ta. 

There exist in the literature several works where this tool has been used in par- 
ticular non-smooth examples, |Hog92| IBK91) . Theorem [3] generalises the result 
stated in [BK91| where the Melnikov method is shown to work, although the 
proof there is not complete. 

The homoclinic version of a piecewise-defined system with a different topology 
was studied in [KunOOj . |Kuk07j and [BF08] , However, the tools developed there 
do not apply for a system of the type (|2.2p . 

We begin by discussing the persistence of objects for e > and r < 1. It 
is clear that by separately extending the systems X^ + eX^ to M 2 x T, where 
T = M./T, we get two smooth systems for which the classical perturbation theory 
holds. It follows then that, as z are hyperbolic hxed points, for e > small 
enough there exist two hyperbolic T-periodic orbits, = {z| t (r);r € [0, T]}, 
with two-dimensional stable and unstable manifolds W s ' u (Af ). 
As the system is non-autonomous, we hx the Poincare section 

6 t0 = {(x,y,t ), (xj)el 2 }, 

and consider the time T stroboscopic map 

nt : 6 to — > @t +T, 

where 

n to (z) = <f>(t + T; t , z, e, r) 

and 4> is defined ^2.41 

This map has zf (to) as hyperbolic fixed points with one dimensional stable and 
unstable manifolds W s ' u (z^(to)) (see Fig. [3]). Proceeding as in |BK91j . we fix 
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the section £ denned in (|2.ip and study its intersection with the stable and 
unstable curves W u (z~(t )) and W s (z^(t j). In the unperturbed conservative 
case (e = and r = 1), W u (z~) and W s (z + ) intersect transversally S in a 
point z . The perturbed manifolds, W u (z~(t )) and W s (z+(t )), intersect T, 
at points z u (to) and z s (to) respectively, e-close to z (see Fig. [5]). Recalling the 
effect of the coefficient of restitution (12.181) explained in i j2.3[ both invariant 
manifolds will intersect if, for some to, one has rz u (to) = z s (to), r < 1. As in 
|BK91j and |Hog92| , we use the unperturbed Hamiltonian H (x,y) to measure 
the distance A(t ,e,r) between z u and z s 

A(t a ,e,r) = H (rz u (t )) - H (z s (t )) = r 2 H (z u (t )) - H (z s (t )). (4.1) 

We then have the following result. 

Theorem 3. Consider system i2.2 \)J2.18\) , and let z = W s (z+) n S. Define 
the Melnikov function 

/+oo 
{H o ,H 1 }{0(t;t o ,z o ,O),t)dt, (4.2) 
-oo 

where 

,(... , n ,_(<f>-(t;t ,zo,Q) ift<t ( .„s 

is the piecewise smooth heteroclinic orbit that exists for r = 1 and e = 0. Assume 
that M(to) possesses a simple zero at to- Then the following holds. 

a) If r = 1, there exists Eq > such that, for every < e < Eq, one can find 
a simple zero = to + 0(e) of the function A(to,e, 1). Hence, the curves 
W u (z~(t ')) and W s {zf{to)) intersect transversally at some point, z^ G S, 
E- close to z a G S aKC? 

{<j)(t;t*,z h ,E),teR}, 
is a heteroclinic orbit between the periodic orbits A~ and Af . 

b) If r < \, there exists p such that, given e, f > satisfying < j < p, one can 
find Sq such that, if £ = eS and r = 1 —rS, then, for < 5 < So, there exists a 
simple zero of the function A(to, sS, 1 — fS) of the form ig = io + 0(j) +0(6) . 
Hence, the curves W u (z~(to)) and W s (zf(to)) intersect E transversally at 
two points, G S, satisfying zt = zo + 0(8) and zZ = zi ' jr, such that 

{^(t;t* ,z+,eS,l-rS),tem} 

is a heteroclinic orbit between the periodic orbits A~ and Af . 

Remark 6. Note that, for r = 1, we recover the classical result given by the 
Melnikov method for heteroclinic orbits extended to the non-smooth system &2.2\) . 
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Proof. Applying the fundamental theorem of calculus to the functions 



s — Ht 



(j> + l-(s;t ,z s,u ,e)), 
we obtain 

H+l- (z s /«) = H±(<f> (r fl /";to,* s/ V) + ±-H+l- (4>+l- (s;t Q ,z s / u ,e) ds) , 

and then make T s l u = +/— oo. However, the limits 

lim </)+/- (t; t , z s / u ,e) 

t— >+/— oo 

do not exist because the flow at the respective stable/unstable manifolds tends 
to the periodic orbit A^. To avoid this limit, we proceed as follows. 
Given to, we define 

U(s) = H- ((f)- {s;t ,z u 7 e)) - H^ ((f)- (s;t , z~ (t ),e)) ,s<t 
f+(s) = H+ (4>+ (s;t Q ,z s ,e)) - Hq (cp+ (s;t Q , z+(t ),e)) , s>t , 

which are well defined smooth functions because the flow is restricted to the 
stable and unstable invariant manifolds or to the hyperbolic periodic orbit and 
never crosses the switching manifold S. 
Then, we write Eq. (|4.1I) as 

A(t Q ,e,r) = r 2 /-(to) - f+(t ) + r 2 H Q - (z~ {to)) - H+ (z+ (t )) . (4.5) 

Noting that 

^(^(to^^+sOT^^U + Ole 2 ), (4.6) 

Cl II 



Eq. (I4.5P becomes 

A(t ,e,r) = r 2 U(t Q )- f + (t ) + (r 2 -l)c 1+ 0(e 2 ) (4.7) 
We apply the fundamental theorem of calculus to the functions (14.41) to compute 

U(t ) = U(T u ) + [° f'_(s)ds = 

f- (T u ) + e jT ( {Ho, fff } (</>' (s; t , z u , e) , s) 
- { H a~, H i} (<t>~ (s\to,z- (t ),e) ,s) jds 



Ufa) = U(T S ) - / f + (s)ds 

Jto 

T" 



(4.8) 



f+(T s )-e / ({H+,H+}(<P+(s;t ,z s ,e),s) 
- {H+, H+} ((f)+ (s; t , z+ (t ) , e) , s) ) ds. 
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Due to the hyperbolicity of the periodic orbits A*, the flow on W S/Au (Ae ) 
converges exponentially to them (forwards or backwards in time). That is, 
there exist positive constants C , A and sq such that 



0+ (s; t , z s ,e)- (j> + (s; t , z+(t ),e) 



< Ce~ Xs , Vs > s , 



and similarly for <fi . This allows one to make T s / U — > +/— oo in Eqs. ()4.8|) . since 

lim /±(s)=0 

S— ¥±00 

and, moreover, the improper integrals converge in the limit. 
Now, expanding the expressions in (|4.8I) in powers of s, we find 

f-(t )=e f° {Ho,Hi}(<f>- ( S ;t ,zo,0),s)ds + O(e 2 ) 

Too (4-9) 

f+(to)=-e / {H+,H+} (<f>+(s;t ,z ,0),s)ds + O{e 2 ), 
Jto 

where we have used property ()4.6p to include the second terms in the integrals 
into the higher order terms. Finally, substituting Eq. (|4.9j) into Eq. (|4.7| . we 
obtain 

A(t , e, r) = (r 2 - l)ci + eM (t ) + 0(e 2 ) + 0(e (r - 1)), (4.10) 

where M(t ) is defined in Eq. P~2"j) . 

We now distinguish between the cases r = 1 and r < 1. If r = 1, we 
recover the classical expression for the distance between the perturbed invariant 
manifolds. By applying the implicit function theorem, it is easy to show that, if 
M(to) has a simple zero at to, then A(to, e, 1) has a simple zero at t^ = to + 0(e). 
Thus, the curves W u (z~ (t^)) and VF s (z+(tg)) intersect E transversally at some 
point, Zh = z w {to) — z s (tj) € S, e-close to z £ S. Therefore, 

{4>(t;t*,z h ,e),t£R}, 

is a heteroclinic orbit between the periodic orbits A~ and A+. 

If r < 1, we define e = e5 and r — 1 — fS, and Eq. (|4.10[) becomes 

Afa ' g y~^ } = -2f Cl + sAf (t„) + 0(<5). (4.11) 

Then we argue as in Theorem [5] As M(to) is a smooth periodic function, it 
possesses at least one local maximum. Let tjvf be the closest value to to where 
M(to) possesses a local maximum, and assume M'(to) ^ for all to between 
to and tu- If M' (to) vanishes between to and tjvf, we then take tu to be the 
closest value to to such that M'(to) = to ensure that M'(to) ^ between to 
and t M - We then define p := ^f^. Then, if 

< ~ < p, 

e 
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there exists to |-close to to such that 

-2fa + M(i ) = 

and M'(to) ^ 0. Hence, we can apply the implicit function theorem to Eq. (|4. 
near the point (to, 8) = (to,0) to conclude that there exists 8q such that, if 
< 8 < Sq, then one can find 

t* =t Q + 0{8) = t Q + O(S) + 0(f/e) 

which is a simple solution of Eq. (|4.11[) . 

Hence, arguing similarly as for r = 1, there exist two points z^ = z s (to) = 
zo + 0(8) and z^ = z u (t^) = z /r + 0{8)r such that z^ = rz^ and 

{<p(t;t* ,z£,eS,l-rS), fel} 

where 

6 (ft* z+ sr) _ir(f,to,t*o,zt/r,e) if t < t* 
Ht ^ Zh ^ r) -\cf>+(t,to,t*,z+,e) ift>*5 
is a heteroclinic orbit between the periodic orbits A~ and A+. □ 

5 Example: the rocking block 

5.1 System equations 

In order to illustrate the results shown in the previous sections, we consider the 
mechanical system shown in Fig. [71 which consists of a rocking block under a 
horizontal periodic forcing given by 

a H (t) = sag cos (fit + 9) . (5.1) 

This system was first studied in |Hou 63 . The equations that govern its be- 
haviour are well known (see for example [YCP80, SK84]). and are given by 

ax + sign(x) sin(a(l — sign (x) x)) = 

— ae cos (a (1 — sign (x) x)) cos (tot) (5-2) 
x(t+) = rx(t- A ) (x = Q) (5.3) 

where the last equation, (|5.3[) . simulates the loss of energy of the block at every 
impact with the ground, as described in £12.31 In addition, the function 



I 1 if ar > 
sign(z) = <^ (5.4) 
1—1 ii x < vi 

distinguishes between the two modes of movement: rocking about the point 
O when the angle is positive (x > 0) or rocking about O 1 when the angle 
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o' o 

a H {t) 



Figure 7: Rocking block 



x is negative. Obviously, this makes the system piecewise smooth and so it 
can be written in the form of Eq. (|2.2j) . Moreover, conditions C.1-C.5 are 
satisfied and, hence, the results shown in previous sections can be applied. 
However, as our purpose here is to illustrate them, we will consider the linearized 
version of Eq. (|5.2p instead, which will permit us to perform explicit analytical 
computations. This linearization is achieved by assuming a -C 1, namely that 
the block is slender |Hog89 . Thus, the system that we will consider, written in 
the form of Eq. ([221), is 



where the perturbation becomes a smooth function due to the linearization. 
If r = 1, system (|5.5[) - (|5.6|) can be written in the form (|2.5|) using the Hamilto- 
nian function 




(5.5) 




(5.6) 



y{t\) = ry{t- A ) (x = o), 



(5.7) 



H E (x,y,t) = H {x,y) +eH 1 (x,t), 



(5.8) 



where 




(5.9) 
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and 

Hi(x,t) = xcoa(u>i) (5.10) 
is T-periodic, with T = 2tt/lu and is a C°° function. 

In addition, when e = 0, conditions C.1-C.5 are fulfilled, and the phase portrait 
for the system (|5.5l) - (|5.6p is equivalent to the one shown in Fig.Q] That is, it pos- 
sesses an invisible fold-fold of centre type at the origin and two saddle points at 
(1, 0) and (— 1, 0) connected by two heteroclinic orbits. Furthermore, the origin 
is surrounded by a continuum of periodic orbits whose periods monotonically 
increase as they approach to the heteroclinic connections. Using Eqs. (|2.15[) 
and (|2.16[) . the symmetries of the Hamiltonian (|5.9[) and assuming yo > 0, these 
periods are given by 



1 



"(2/0 ) = 4 / —j=====dx 

Jo y/yfi + x z - 2x 



= 21nl'i±^L (5.11) 



i - yo J ' 



and hence a'(y ) > 0. 



5.2 Existence of periodic orbits 

We first study the persistence of (n, m)-periodic orbits for r = 1 in Eq. (|5.7|) by 
applying Theorem [TJ The subharmonic Melnikov function (|3.3[) becomes 

pnT 

M n > m (t ) = - / U y (q c (t)) cos(w (f + t ))dt, (5.12) 



where q c (t) is the periodic orbit of the unperturbed version of system f|5.5[) - ()5.6[) 

-2 

with Hamiltonian c = -y- satisfying q c (0) = (0,y~o) and 

, fnT\ eS - 1 

?7o = a" 1 — = — • (5.13) 

\m J e ^ + 1 

We now want to obtain an explicit expression for Eq. (|5.12j) . Thus we first note 
that the solution of system (|5.5[) - (|5.6[) with initial condition (xq, yo) at i = to is 
given by 

x ± {t) = Cfe t + Cfe~ t ±l (5.14) 
y ± (i) = C 1 ± e t -C* 2 ± e- t , (5.15) 



where 



As explained in £12.41 the superscript + is applied if xq > or xq — and yo > 0, 
and the — otherwise. 



Cf = X0+ l° T V *°, Ct = X °~^ 1 e*°. (5.16) 
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Assuming xq — and yo = yo > 0, an expression for H y (q c (t)) becomes 

{nT 
C 1 e t -C 2 e-\ if0<t< — 
nT „T < 5 - 17 ) 

- Cl e*-^ + C 2 e-^, if^<i<^, 
2m to 

where 

Cl = ^,C 2 = ^f±. (5.18) 
Thus, Eq. (pTT2"j) becomes 

M"- m (t ) = -|] (C ie 4 - C 2 e- 4 ) cos ^ (t + to+j^yjdt 

+ ^-Cie*"^ + C 2 e~ t+ ^ cos (* + *o+j^- 



and, after some computations, we have 

4 



cos (cuto) , if m = 1 



M n ' m (t )={ Lo 2 + l'^ y "" J " - (5.19) 

0, if to > 1. 

As M n,1 (io) has two simple zeros, to = -j and ^ = by Theorem [U if 
e > is small enough, the non-autonomous system (|5.5l) - (|5.6p possesses two 
subharmonic (n, l)-periodic orbits. In addition, the initial conditions of these 
periodic orbits are e-close to 

(0,y ,t l ) = (Q,^—L-) (5.20) 
e 2 + 1 4 

and 

e^f" — 1 3T 

(0,y Ol ^) = (0,^— ,— ), (5.21) 
e 2 + 1 4 

respectively. 

Proceeding as in Remark [3l one can solve numerically Eq. (|3.2[) with to = 1 
and find the initial conditions for such a periodic orbit. In Fig. [8] we show the 
result of that for n = 5. Each periodic orbit is obtained by using the points 
given in Eqs. (15.20)) and (|5.21[) to initiate the Newton method. Then, tracking 
the obtained solution, e has been increased up to e = 1.6565 • 10~ 2 . 

Regarding the existence of (n, m)-periodic orbits with m > 1 (ultrasubhar- 
monic orbits), as the subharmonic Mclnikov function is identically zero nothing 
can be said using the first order analysis that we have shown in this work. 
However, if instead (|5. 10[) one considers the perturbation 

H\(x,t) = x (cos (wt) + cos (kujt) ) , 
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Figure 8: Periodic orbits for n — 5 and m = 1, oj = 5 and e = 1.6565 • 10 2 . 
Their initial conditions are e-close to the points given in Eqs. f|5 . 20[) and (|5.2ip . 



then, it can be seen that the corresponding Melnikov function possesses simple 
zeros for m = k and n relatively prime odd integers. Thus, periodic orbits im- 
pacting m > 1 times with the switching manifold can exist if higher harmonics 
of the perturbation are considered. 

Let us now introduce the energy dissipation described in £13.21 and con- 
sider the whole system (|5.5|) - (|5.7j) with r < 1 using the Hamiltonian pertur- 
bation (|5.8j) . From Theorem [21 simple zeros of the Melnikov function (|5.19[) 
also guarantee the existence of (n, l)-periodic orbits when 1 — r is small enough 
with respect to e. More precisely, taking 

£ = id, r = l-?8, (5.22) 

condition (|3.14[) becomes 

< ~- < \ ( e l T +1 ) M n >\t M ) := p, (5.23) 

e 2 y e — - i j 

where M ra,1 (tjvtf-) = M n,1 (-j) = u * +i is the maximum value of the Melnikov 
function (|5.19p . Then there exists an (n, l)-periodic orbit if 5 > is small 
enough. The initial condition of the periodic orbit is located in a ^-neighbourhood 
of the point (a^J/cb^o) = (0,2/Oj*o)j where yo is defined in Eq. (I5.13|) . such that 

a(£o) = nT 

and to is given by the simple zeros of Eq. (j3. 13[) , which becomes 

-2fyl + iM n ' 1 {t ) = Q. (5.24) 
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(a) S = 1.04 



(b) S = 4.125 



Figure 9: (5, l)-periodic orbits for cj = 5 and Z = 0.07. Tracking the obtained 
solution, the perturbation parameter S has been increased up to its maximum 
value. Initial conditions close to (yo,^) and (j/oi^o) have been used in (a) and 
(b), respectively. 



Hence we find 
3 1 




$ = ~arccos — - +(<-!)_, < = 1, 2. (5.25) 



As before, we set n = 5 and to = 5. Then expression (|5.23p becomes 

< - < 0.0914. (5.26) 

£ 

Hence, for any fixed ratio | satisfying (|5.26|) there exist two points, (j/o,io), 
i = 1,2, such that, if S is small enough, Eq. (|3.10j) possesses a solution 5-close 
to them. Such a solution is an initial condition for an (n, l)-periodic orbit of 
system (|5.5[) - (I5.7I) . with r = 1 — rd and e — ES. 

In Fig. [9] some of these orbits are shown for one value of the ratio satisfy- 
ing (|5.26p . Two different periodic orbits are shown, which are the ones whose 
initial conditions are <5-close to the values (yo,to) and (yo,t^)- In both cases, S 
tracks the solution, up to the value from which solutions of Eq. (|5.23|) can no 
longer be found. These are the values used in the simulations shown in Fig. |H1 
Above the limiting value of the ratio given in (15.26[) . no (5, l)-periodic orbits 
were found for uj = 5. 



5.3 Existence curves 

We now use Theorem[2]to derive existence curves for the (n, l)-periodic orbits (n 
odd) and compare them with the ones obtained in (Hog89 . Unlike in [Hog89 , 
we obtain these curves in the r-e plane through the variation of S. 
The limiting condition provided by Theorem [2] is given in Eq. (|5 .23[) . Thus, for a 
given r close to 1 (rS close to 0), it is natural to arbitrarily fix r and minimize e 
by maximizing the ratio in (|5.23p . setting ? = p. However, the upper boundary 
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of 6, So, provided by Theorem [5] tends to zero as | — > p, as it is derived from 
the implicit function theorem. Thus, it is not possible to uniformally bound S 
for all the ratios between and p. Hence, the condition | = p can not be used 
to derive a limiting relation between r and e. Instead, we proceed as follows. 
We first fix n odd and u > 0. Then, for every ratio < | < p, we increase 5 
from to 5q by numerically tracking the obtained solution using as initial seed 
one of the values provided in Eq. (|5.20|) or (15.211) . This results in a curve in the 
r-e plane parametrized by the rati o Z. 

Regarding the results obtained in |Hog89| , such a curve was obtained analyti- 
cally and for global conditions, and has the expression 

, , (l+ W 2 )i?(l-cosh(^)) , 

emin(fl) - , V ' V V2JJ = , (5.27) 

^lj 2 sinh 2 i? 2 + (2 - R) 2 (1 + cosh (^)) 2 

where R = 1 — r. 

As our result is only locally valid, in order to compare both results we have to 
check whether both curves are tangent at e = 0. From (|5.27|) we easily obtain 

2 \e— + IJ P 
which, by the inverse function theorem, tells us that both curves are tangent at 

£ = 0. 



In Fig. [TUl we show an example for n = 5 and to = 5 using initial conditions 
near (I5.20[) . As can be seen, the curve provides, for every value of r, both the 
maximum and minimum values of e for which a (5, l)-periodic orbit exists. The 
lower boundary derived in |Hog89| , (e m in( - )) ^ 1 ( £ ) i s a ^ so shown. As demon- 
strated above, both curves are tangent at e = 0, with slope equal to p. Note 
that the lower boundary does not coincide with the line 1 — r = pe, although 
their difference tends to zero as r — > 1. This confirms that one can not derive 
the minimum value of e from condition (|5 . 23[) for every fixed r. 
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